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Results on the production of the double-strange cascade hyperon H” are reported for collisions of 
p(3.5 GeV) + Nb, studied with the High Acceptance Di-Electron Spectrometer (HADES) at SIS 18 at GSI 
Helmholtzzentrum for Heavy-Ion Research, Darmstadt. For the first time, subthreshold production is ob¬ 
served in proton-nucleus interactions. Assuming a phase-space distribution similar to that of A hyperons, 
the production probability amounts to Pg- = (2.0 ±0.4 (stat) ±0.3 (norm) ±0.6 (syst)) x 10“* resulting in 
a H“/(A ± S®) ratio of Pg- / Pa+eo = (1-2 ± 0.3 (stat) ± 0.4 (syst)) x 10“®. Available model predictions 
are significantly lower than the estimated yield. 

PACS numbers: 25.75.Dw, 25.75.Gz 

The double-strange S' baryon (also known as cascade par- strangeness conservation, NN —> NSKK. In fixed-target ex- 

ticle) when produced in elementary nucleon-nucleon (NN) periments, this requires a minimum beam energy of Ethr = 

collisions must be co-produced with two kaons ensuring 3.74 GeV (y^Sthr = 3.25 GeV). In heavy-ion and even in 
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nucleon-nucleus collisions cooperative processes are possi¬ 
ble allowing for the production below this threshold. Above 
threshold and in heavy-ion reactions, the S” hyperons were 
measured over about three orders of magnitude of the centre- 
of-mass energy covered by the LHC (s/snn = 2.76 TeV 
El), RHIC = 62.4, 200 GeV HHEl), SPS = 

8.9,17.3 GeV El, = 6.4 - 17.3 GeV) S), and 

AGS {y/sNN = 3.84 GeV lO) accelerators. The yield of 
multi-strange particles produced in nucleon-nucleus (p-tA) 
and nucleus-nucleus (A + A) collisions below their produc¬ 
tion threshold in NN collisions, is expected to be sensitive 
to the equation of state (EoS) of nuclear matter, similar to 
single-strange hadrons GHSl. In heavy-ion reactions, the nec¬ 
essary energy for the production of multi-strange hyperons 
can be accumulated via multiple collisions involving nucle¬ 
ons, produced particles and short-living resonances. The cor¬ 
responding number of such collisions increases with the den¬ 
sity within the reaction zone the maximum of which in turn 
depends on the stiffness of the EoS. 

Predictions of sub-threshold cascade production at ener¬ 
gies available with the heavy-ion synchrotron SIS 18 at GSI, 
Darmstadt, were made within a relativistic transport model 
EHl. The cross sections of the strangeness exchange reac¬ 
tions KY —ttS (Y = A, E), which were thought to be es¬ 
sential for 5 creation below NN threshold, were taken from 
a coupled-channel approach based on a flavor SU(3)-invariant 
hadronic Lagrangian m. At that time, no subthreshold 5 
production was observed; the first announcement came from 
the HADES collaboration studying, at SIS 18, Ar-tKCl re¬ 
actions at a beam kinetic energy of 1.76A GeV (y/s^N = 
2.61 GeV ESI)- The deduced 5“/A ratio was found substan¬ 
tially larger than any model prediction available at that time. 
Shortly after, other strangeness-exchange reactions, e.g. the 
hyperon-hyperon scattering processes, YY —5N, exhibiting 
quite high cross sections, were figured out to largely account 
for this discrepancy E3l [T4ll . while the reaction KY —ttS 
was found negligible ES- Also, a very recent investigation 
ES of deep-subthreshold 5 production in nuclear collisions 
with the UrQMD transport model ESEl making use of the 
YY cross sections provided in ref. E3l showed that the hy- 
peron strangeness exchange is the dominant process contribut¬ 
ing to the S yield. However, the model could not satisfactorily 
explain the S“/A ratio measured by HADES. Presently, no 
experimental data exist on S production in p H- A interactions 
near threshold. The lowest energy so far at which 5“ produc¬ 
tion has been observed in collisions of p H- Be and p H- Pb is the 
maximum SPS energy (y/s^N = 17-3 GeV ESI)- 

It would be interesting to learn which processes contribute 
mainly to subthreshold S production in case of nucleon- 
nucleus collisions, which are considered as a link between 
elementary NN and heavy-ion collisions. Eor instance, in 
pH- A reactions at a beam kinetic energy of 3.5 GeV (fixed 
target, y/sNN = 3.18 GeV), the scattering of two incoher¬ 
ently produced hyperons appears rather improbable. Also, at 
first sight, direct double-hyperon production seems to be im¬ 
possible, since the threshold of y/Sthr,AA = 3.22 GeV for the 


channel requiring the lowest energy effort, pp —> AAK+K+, 
is only marginally lower than the 5 threshold. However, al¬ 
ready the consideration of a rather modest Eermi motion of 
the nucleons within the nucleus, i.e. a counter-motion against 
the projectile with a momentum of about 50 MeV/c (being 
well below the Eermi momentum), would lift the available 
energy above both thresholds. Another possibility to gain 
the necessary energy would be the scattering of the proton 
at an object acting more massively than a single nucleon in 
the nucleus. Thus, for p (3.5 GeV) + A collisions, a target ob¬ 
ject X as heavy as only 1.11 nucleon masses is sufficient to 
reach the threshold energy of the final state X5KK. Hence, 
the co-operation of two, or more, correlated target nucleons, 
e.g. bound in a particles, would allow for this kinematic ef¬ 
fect. Correspondingly, a high collectivity of the target nucleus 
was already observed in deep-subthreshold kaon production 
in p H- A (A = C, Cu, Au) collisions at a beam kinetic energy 
of 1.0 GeV studied by ANKE at COSY-Jiilich EH. Also, in 
electron scattering experiments, ^^C(e,e’p) at 4.627 GeV, per¬ 
formed at JLab ll20l . a surprisingly high fraction of nucleons, 
i.e. 20 %, are found to be strongly correlated, predominantly 
in the form of proton-neutron pairs. Einally, 5 production 
in p H- A (and A H- A) might happen via an exotic channel as 
pA —> 5K, probably via the excitation and decay of a massive 
resonance as, e.g., A*(2000, 2100,...) IJT]. Another possible 
scenario is to produce a non-strange heavy resonance with a 
mass sufficiently high to decay into SKK. This hypothesis is 
motivated by the fact that in pH-p collisions at 3.5 GeV beam 
kinetic energy a substantial fraction of the exclusive produc¬ 
tion of S(1385)+-b K+ -b n was found to proceed via an in¬ 
termediate broad A++ excitation at about 2000MeV fT2\ . 

In this Letter, we report on the first observation of sub¬ 
threshold 5” production in nucleon-nucleus collisions at 
y/sNN — y/sthr = — 70MeV. The experiment was per¬ 
formed with the High Acceptance Di-Electron Spectrometer 
(HADES) at the Schwerionensynchrotron SIS 18 at GSI, 
Darmstadt. HADES, although primarily optimized to measure 
di-electrons ll23l . offers also excellent hadron identification 
capabilities Il24lj26l . A detailed description of the spectrome¬ 
ter is presented in ref. lEil. The present results are based on a 
dataset which was previously analyzed with respect to e+e~ 
ll28ll as well as to pion and 77 1291 . K° 1^ and A 1311 produc¬ 
tion in collisions of pn-Nb at 3.5 GeV. The main features of 
the apparatus relevant for the present analysis are summarized 
in ref. l3Ti . 

In the present experiment, a proton beam of about 2 x 10 ® 
particles per second with kinetic energy of 3.5 GeV was inci¬ 
dent on a 12-fold segmented target of natural niobium (®®Nb). 
The data readout was started by different trigger decisions 
I29I . Eor the present analysis, we employ only the data of the 
first-level (LVLl) trigger, requiring a charged-particle multi¬ 
plicity > 3 in the time-of-flight wall composed of plastic scin¬ 
tillation detectors. We processed about Alvli = 3.2 x 10 ® of 
such LVLl events. 

It is important to mention that E® hyperons decay almost 
exclusively into A’s via the decay E® —> A 7 (branching ratio 
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BR = 100%, CT = 2.22 X 10“^^ m Erl '), with the photon 
not being detected in the present experiment. Hence, through¬ 
out the paper, any “A yield” has to be understood as that of 
A + YP. Correspondingly, in case of simulations, where the 
individual particle species are known, the yields of A and 
hyperons are summed up. 

In the present analysis, we identify the and A hyperons 
through their weak decays > Att” {BR = 99.9 %, cr = 
4.91cm) and A —> p7r“ {BR = 63.9%, ct = 7.89 cm) ||2TI . 
with the charged hadrons detected in HADES lfT2l |2^ [311 . 
The long lifetimes cause a sizeable fraction of these particles 
to decay away from the primary vertex. The precision of the 
track reconstruction with HADES is sufficient to resolve these 
secondary vertices na. To allow for A selection various topo¬ 
logical cuts on single-particle and two-particle quantities were 
applied. These are i) a minimum value of the proton track^ 
distance to the primary vertex (p-VecToPrimVer), ii) the same 
for the 7r“ (7ri-VecToPrimVer), iii) an upper limit of the p-7r“ 
minimum track distance (p-TTi-MinVecDist), and iv) a min¬ 
imum value of the A decay vertex distance to the primary 
vertex (A-VerToPrimVer). Here, the off-vertex cut iv) is the 
main condition responsible for the extraction of a A signal. 
Starting with the moderate conditions as used in the previous 
high-statistics analysis of the A phase-space distribution and 
polarization (ED, a clear A signal could be separated from 
the combinatorial background in the p-7r” invariant-mass dis¬ 
tribution. While in that analysis a signal-to-background ra¬ 
tio in the order of unity was sufficient, for the present 5” 
search we start with a higher A purity (>85 %, cp. flEil l. 
Hence, with the stronger cuts and the requirement of an ad¬ 
ditional 7r“ meson, the number of reconstructed A hyperons 
decreases from about 1.1 million to 300,000. (No event con¬ 
taining clearly more than one A was found.) Taking this still 
high-statistics A sample, we started the investigation by 
combining - for each event containing a A candidate (selected 
by a ±2 ct window around the A peak) - the A with those tt~ 
mesons not already contributing to the A. The result was a 
structureless A-tt” invariant mass distribution. Hence, ad¬ 
ditional conditions were necessary: v) a lower limit on the 
2 nd 7r“ (potential 5“ daughter) track distance to the pri¬ 
mary vertex (7r2-VecToPrimVer), vi) an upper limit of the 
distance of the pointing vector w.r.t. the primary vertex 
(5-VecToPrimVer), vii) a maximum value of the minimum 
track distance of the A and the 2nd 7r“ (7r2-A-MinVecDist), 
and viii) a minimum value of the distance of the vertex 
relative to the primary one (5-VerToPrimVer). 

Starting with the cut settings used in our previous analysis 
of deep-subthreshold production in collisions of Ar + KCl 
at 1.76A GeV 02 and optimizing further for the present 
experiment which exhibits different multiplicities and phase- 
space distributions of the involved particles, we find a signifi¬ 


' With “track” we mean the trajectory of a particle track extrapolated up to 
the relevant vertex. 


cant narrow peak structure in the A-tt” invariant-mass distri¬ 
bution displayed in Eig.[2 (Eor convenience, we use identical 
mass and energy units.) The position is slightly lower by about 
4 MeV than the PDG value of 1321.71 MeV ETIl . very prob¬ 
ably due to a minor systematics uncertainty of the momentum 
calibration for charged particles in the inhomogeneous field 
of the toroidal magnet which leads to a slight (<0.4 %) phase- 
space dependence of the mass of the reconstructed weakly 
decaying mother particle as observed already for K° and A 
reconstruction IMIHI- The width of the present peak, how¬ 
ever, is well in agreement with the results of GEANT ll32]l 
simulations of about 2-3 MeV for A and 5” hyperons. More 
importantly, also the cut dependences of the yield of the po¬ 
tential 5“ match well those found from GEANT simulations 
(see below). Thus, we attribute the signal to the decay of the 
5“ baryon. 

The full curve in Eig.[2 shows the result of a fit to the data 
with a model function consisting of a Gaussian function for 
the peak and a polynomial function of 2nd order for the com¬ 
binatorial background (bg). Integration around the peak max¬ 
imum within a window of ±5 MeV (« ±2(T, with cr being 
the Gaussian width) we find N-^- = 90 ± 18 with the sta¬ 
tistical error given. The signal-to-background ratio and the 
significance, signal/-\/signal + bg, amount to 0.39 and 5.0, 
respectively. Note that the raw 5“ yield per LVLl event of 
-^H-/-^LVLi = 2.8 X 10“® is yet a factor seven smaller than 
the corresponding yield in Arn-KCl reactions at 1.76A GeV 
d. We studied also the raw phase-space distribution of the 



[MeV] 


FIG. 1: The experimental A — 7r“ invariant-mass distribution. The 
error bars show the statistical errors. The curve represents a combi¬ 
nation of a Gaussian and a polynomial function used to fit the data. 


baryons. To that purpose, the yield within a window of 
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±5 MeV around the S“ peak in Fig.was selected, and the 
combinatorial background below the peak was subtracted with 
the help of a corresponding sideband analysis. The result¬ 
ing transverse-momentum vs. rapidity distribution was found 
strongly biased by the HADES acceptance, i.e. essentially by 
the lower and upper polar angle limits of 18 and 85 degrees 
(2911311. This finding is confirmed by studies of the detector 
acceptance of simulated data and found to be rather indepen¬ 
dent of the input phase-space distributions. The mean value 
and the r.m.s. width of the experimental rapidity distribution 
amount to 0.54 and 0.16, respectively. The corresponding val¬ 
ues of the transverse-momentum distribution are 0.52 GeV/c 
and 0.17 GeV/c, respectively. 

Corrections for the finite acceptance and reconstruction ef¬ 
ficiency were deduced from simulations. Thermo-statistically 
distributed 5“ baryons, characterized by a temperature pa¬ 
rameter T, were generated with the event generator Pluto 
Since the phase-space distribution of the is not 
known, the experimental A phase-space distribution (found 
to be strongly influenced by hyperon-nucleon collisions (31]) 
served as benchmark for the S” hyperon. Consequently, in 
Pluto we allowed for longitudinally shifted and elongated 5” 
phase-space distributions. For this purpose, two longitudinal 
shape parameter, i.e. the mean, (y) = 0.3, and the width, 
ay = 0.57, following from a Gaussian fit to the A rapid¬ 
ity distribution are introduced. We investigated the 
geometrical acceptance for a broad range of transverse and 
longitudinal shape parameters, i.e. T = 50,65,80, 95 MeV 
(cf. ref. Oil), (y) = 0,0.3,0.6. With the given parame¬ 
ters, we determined, with Pluto, the average HADES accep¬ 
tance for the 5” hyperon (including the branching ratio of 
64% for the decay of its daughter, A —p7r“) and its varia¬ 
tion within the parameter ranges. Thus, we estimated a sys¬ 
tematic error of about ±25 % around the average, purely ge¬ 
ometric, acceptance of Cacc, sym = 6.4 • 10“^ for the 
above given phase-space parameters. The same Pluto data are 
processed through GEANT, modeling the detector response. 
The GEANT data were embedded into real experimental data 
and then processed through the full analysis chain (using the 
same topological cuts i) ... viii) as applied to the experimen¬ 
tal data). The mean value and the r.m.s. width of the result¬ 
ing HADES-filtered 5” rapidity distribution amount to 0.60 
and 0.16, respectively, quite similar to the experimental val¬ 
ues. Relating the output to the corresponding input, the total 
acceptance x reconstruction efficiency was estimated to 
eggf = (8.49 ± 0.24) • 10“®. As in the A hyperon analysis 
(111 we correct for the LVLl trigger bias w.r.t. minimum- 
bias events, Flvli = A^minbias/A^LVLi = 1-53 ± 0.02, and 
for empty-track events due to non-target interactions, Fmt = 
0.17. Finally, we note that the experimental A rapidity distri¬ 
bution on does not appear perfectly symmetric. Also, trans¬ 
port model calculations (cf. Fig. 6 of ref. On i rather predict 
a faster yield decrease in the backward hemisphere than is ex¬ 
pected from the rapidity-symmetric Pluto distribution. Pro¬ 
vided that the hyperon exhibits a similar asymmetric ra¬ 
pidity distribution as the A, we have to correct the yield 


for this difference. We do that by calculating the acceptance 
ratio of an asymmetric rapidity distribution to the symmetric 
one. The asymmetric distribution is modelled by a function 
consisting of two Gaussian distributions. The first Gaussian is 
the above one describing the A rapidity density in the region 
where experimental data points are available and the second 
Gaussian is a more narrow one (ay = 0.13, (y) = 0.07) fol¬ 
lowing the transport model predictions at lower rapidities and 
joining up to the first Gaussian at y ~ 0.06. The ratio of the 
acceptances for asymmetric and symmetric rapidity distribu¬ 
tions amounts to Easy ~ t^acc, asy/t^acc, sym “ 1.32 ± 0.02. 
Assuming that the total acceptance and reconstruction ef¬ 
ficiency can be factorized into an acceptance part and and a 
pure reconstruction-efficiency part which itself does not vary 
within the, rather limited, acceptance, this factor is used to 
correct Ceff. 

With all the necessary quantities and correction factors at 
hand, we calculate the production probability to 

P,_=._ ^ __ 

(1 — Amt) Alvli A^lvli Aasy CeS 

= (2.0 ± 0.4 (stat) ± 0.3 (norm) ± 0.6 (syst)) x 10“"^, 

( 1 ) 


where the statistical, absolute normalization, and systematic 
errors are given. 

The dependences of the yield on the cut values of the 
various geometrical quantities are displayed in Fig.[^ Due to 
the limited 5“ statistics, only one cut could be varied while all 
the others are kept fixed to the optimum values (indicated by 
arrows) yielding the most significant signal. The dependences 
of experimental data (full circles) and GEANT simulations 
(open circles) are found to be in good agreement. 

Taking the A production probability of 0.017 per minimum- 
bias event as estimated in ref. (32, the ratio of ^ and A pro¬ 
duction yields can be determined. Such a ratio, when derived 
from the same data analysis, has the advantage that system¬ 
atic errors (e.g. the uncertainty of the absolute normalization) 
cancel to some extent. The ratio amounts to 


= (1.2 ± 0.3 (stat) ± 0.4 (syst)) x lO'^. (2) 

”a+so 


Here, the statistical error is dominated by the 20 % error of 
the signal, while the systematic error is governed by the 
stability of the signal against cut variations and by the range 
of the parameters entering the simulation. 

The deduced ratio Q can be compared with correspond¬ 
ing ratios at higher energies lfTI46l ITSl l34l - l3^ . Eigurej^ shows 
a compilation of such ratios as a function of ^ysN'N■ So 
far, the lowest energy at which a S“/A ratio is available is 
■y/SNN = 2.61 GeV, i.e. 630 MeV below the threshold in 
NN collisions. The corresponding ratio (open circle) was ex¬ 
tracted by HADES from Ar + KCl reactions at a beam energy 
of 1.76A GeV ifT^ . A steep decline of the 5“/A produc¬ 
tion ratio is observed around threshold, where now a sec¬ 
ond data point (full circle) is available at an excess energy 
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Tij-A-MinVecDist E-VerToPrimVer 


FIG. 2: Relative yield as a function of the cut value of various A 
and geometrical distances (see text, abscissa units are mm). The 
full (open) circles display the experimental (simulation) data. The 
vertical and horizontal arrows indicate the chosen cut values and the 
region of accepted distances, respectively. 


of -70 MeV. This allows for comparisons to model calcula¬ 
tions (see below). To visualize the energy dependence of the 
proton-induced data (full curve in Fig.[^, we fitted the corre¬ 
sponding ratios with a function/(x) = C{1—{D/x)^Y {W\\h 
X = C = 0.44, D = 2.2 GeV, /r = 0.027, v = 0.78), 

a simple parameterization which may be used to estimate the 
expected /A ratio in energy regions, where data are not yet 
available. 

The S“/(A -F S°) ratio has been investigated within a sta¬ 
tistical approach. We performed a calculation with the pack¬ 
age THERMUS Ell, using the mixed-canonical ensemble, 
where strangeness is exactly conserved, while all other quan¬ 
tum numbers are conserved only on average by chemical po¬ 
tentials. The optimum input parameters for this calculation 
(i.e. temperature, T = (121 ± 3) MeV, baryon chemical po¬ 
tential, /rs = (722 ± 85) MeV, charge chemical potential, 
PQ = (24 ± 20) MeV, fireball radius, R = (1.05 ± 0.15) fm, 
and radius of strangeness-conserving canonical volume, Rc = 
(0.8±2.1) fm) follow from the best fit to the available HADES 
particle yields (7r“, 7r°, p, w, K°, A) in p + Nb collisions at 
3.5 GeV II 29 H 3 II . We obtained a 5" yield of 1.0 x 10"^ and 
a S“/(A -F E°) ratio of 8.1 x 10“"* (asterisk in Eig.[^. Both 
values are significantly lower than the corresponding experi¬ 
mental data. 

We also estimated the 5 production probability within two 
different transport approaches, both having implemented the 
aforementioned strangeness-exchange channels. The first ap- 



FIG. 3: The yield ratio H“/(A -F S'^) as a function of ^^snn or 
\/snn — \/sthr (inset). The arrows indicate the threshold in free NN 
collisions. The open symbols represent data for symmetric heavy-ion 
collisions measured at LHC fni34l (cross), RHIC |2l|3 (stars), SPS 
Em (triangles), AGS O (square), and SIS 18 1121 (circle). The 
filled cross depicts p + p collisions at LHC on, while the downward 
and upward pointing filled triangles are for p + A reactions at DESY 
(36) and SPS EH, respectively. The filled circle shows the present 
ratio for p (3.5 GeV) -l-Nb reactions (statistical error within ticks, 
systematic error as bar). The full curve is a parameterization (see 
text) of the proton-induced reaction data. The asterisk, diamond 
and filled star display the predictions of the statistical-model pack¬ 
age THERMUS 03, the GiBUU OHISIl, and the UrQMD fTtl[T3 
transport approaches, respectively. 


proach is the UrQMD model Id 113 (version^ 3.4). For 5 
hyperons, we derived a yield of (6.9 =F 2.8) x 10“^ per event 
which is more than two orders of magnitude lower than the 
experimental yield 0 and decreases only by a factor of two, 
if the channels YY —SN (with cross sections from ini) 
are deactivated; i.e. in the model hyperon-hyperon fusion 
is of minor importance for S production in proton-nucleus 
reactions at 3.5 GeV. The A rapidity distribution, however, 
was fairly well reproduced by UrQMD BTI . The resulting 
S“/(A-FS°) ratio amounts to (3.1 =F 1.2) x 10“® (filled 
star in Fig.[^. The second transport approach we used is the 
GiBUU model ll38l [39l (release^ 1.6). We estimated a 5” 
yield of (6.2 zF 0.9) x 10“®, a value being considerably higher 
than the prediction by the UrQMD model, but still signifi¬ 
cantly lower than the experimental yield 0 - Also here, the 


^ http://urqmd.org 
^ https://gibuu.hepforge.org 
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total A yield was quite well (up to 90 %) reproduced OTl . The 
S“/(A + S°) ratio amounts to (3.8 ±0.5) x 10“^ (filled dia¬ 
mond in Fig.0. The difference of the results of both transport 
models may origin from different parameterizations of cross 
sections of elementary processes. 

Summarizing, we investigated the production of the hy- 
peron in collisions of p(3.5 GeV)H-Nb. For the first time, 
subthreshold S production is observed in proton-nucleus in¬ 
teractions. Assuming a phase-space distribution similar 
to that of A hyperons, the yield per event amounts to 
(2.0 ±0.4(stat) ±0.3 (norm) ±0.6(syst)) x 10“"*. Tak¬ 
ing advantage of a recent investigation of A hyperon produc¬ 
tion and polarization in the same collisions system ll3Tl . the 
S“/(A ± S°) yield ratio of (1.2 ± 0.3 (stat) ± 0.4 (syst)) x 
10“^ is derived. Corresponding estimates with the statistical- 
model package TFIERMUS are significantly lower than the 
experimental data, by more than an order of magnitude. The 
GiBUU transport approach predicts a yield of similar 
level as that of THERMUS. The UrQMD transport model ex¬ 
tremely underestimates the present data, i.e. the yield is an 
order of magnitude lower than the GiBUU result. Both trans¬ 
port codes, however, reproduce the pion and single-strange 
hadron yields fairly well. Hyperon-hyperon scattering pro¬ 
cesses, YY —> SN, recently accounted for the temporary puz¬ 
zle of deep-subthreshold S production in heavy-ion reactions, 
are found to be of minor importance for subthreshold S gen¬ 
eration in proton-nucleus collisions. Hence, a new S puzzle 
appears, now in proton-nucleus collisions. Probably, direct 
subthreshold S production happens on the tail of the Eermi 
distribution of the nucleons within the nucleus, or some of the 
target protons and neutrons form strongly correlated nucleon 
pairs, kinematically allowing for either direct S production 
or the population of an intermediate massive resonance de¬ 
caying into S baryons. Consequently, more systematic inves¬ 
tigations are necessary, including nucleon-nucleon and pion- 
induced reactions, goals already pursued by HADES. 
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